Microscopic dynamics and relaxation processes in liquid Hydrogen Fluoride 
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Inelastic x-ray scattering and Brillouin light scattering measurements of the dynamic structure 
factor of hquid hydrogen fluoride have been performed in the temperature rangeT = 214 2S?iK. 
The data, analysed using a viscoelastic model with a two timescale memory function, show a positive 
dispersion of the sound velocity c(Q) between the low frequency value co(Q) and the high frequency 
value Cooa(Q). This finding confirms the existence of a structural (a) relaxation directly related to 
the dynamical organization of the hydrogen bonds network of the system. The activation energy 
Ea of the process has been extracted by the analysis of the temperature behavior of the relaxation 
time Ta{T) that follows an Arrhenius law. The obtained value for Ea, when compared with that 
observed in another hydrogen bond liquid as water, suggests that the main parameter governing the 
Q-relaxation process is the number of the hydrogen bonds per molecule. 

PACS numbers: 61.20.-p, 63.50.+X, 61.10.Eq, 78.70.Ck 



I. INTRODUCTION 



To understand how the presence of a relaxation pro- 
cess affects the dynamics of the density fluctuations in 
liquids is one of the open problems in the physics of the 
condensed matter. Despite the fact that a large effort 
has been devoted to shed light on this subject, the mat- 
ter is still under debate. In this respect, among all the 
relaxations active in a liquid, particular attention has 
been paid to relaxation processes of viscous nature which 
strongly affect the longitudinal density modes. They in- 
clude at least two distinct contributions: a structural (or 
a) and a microscopic (or /i) process. The a-process is as- 
sociated to the structural rearrangement of the particles 
in the liquid and its characteristic time (tq) is strongly 
temperature dependent, can vary several order of 
magnitude going from the ps, in the high temperature 
liquid phase, to ~ 100 s in glass-forming materials at 
the glass transition temperature. The /i-process takes 
its origin from the oscillatory motion of a particle in the 
cage of its nearest neighbors before escaping. Its char- 
acteristic time (r^) is shorter than and its "strength" 
is often larger than the strength of the a-process. Other 
relaxation processes, beyond the a and the instantaneous 
processes, associated with the internal molecular degrees 
of freedom may be observed in molecular liquidsS^iSi. 
The existence of the a and /i processes, already intro- 
duced several years ago in a molecular dynamic simula- 
tion study on a Lennard- Jones fluidi, has recently been 
proved by experiments on liquid metalsSiiii^. In this 
respect another very important class of liquids to con- 
sider are the hydrogen bonded (HB) liquid systems. In 
these compounds indeed, the highly directional hydrogen 
bond plays a crucial role in the determination of their 
microscopic properties. Despite the large number of the- 



oretical studies 6/18^9^10, 11^12,13, 14^ -(.j^g ^g^y which the 
peculiarities of the hydrogen bond networks affect the 
static organization and the dynamical behavior of these 
compounds is still a subject of discussion. Many are the 
parameters related to the hydrogen bond that must be 
considered in the description of the physical properties of 
these liquids, as for example the hydrogen bond strength, 
the spatial network arrangement of the hydrogen bonds 
and the number of hydrogen bonds per molecule. From 
an experimental point of view, a study of the collective 
dynamics as a function of these parameters is extremely 
important to clarify the role played by each of them on 
the physical properties of HB systems. Among the HB 
liquids, hydrogen fluoride (HF) represents one of the most 
intriguing systems as demonstrated by the large amount 
of theoretical study on its static ^iSiSiiSiii and dynamic 
-'^^i^'^i-'^'^ properties. It represents, in fact, a perfect HB 
model system: it has a simple diatomic molecule and a 
very strong hydrogen bond that determines a linear chain 
arrangement of the HB network. Nevertheless, despite its 
apparent simplicity, only few experimental data are avail- 
able because of the very high reactivity of the material 
that consequently makes its handling extremely difflcult. 
In a previous work^^ we studied the high frequency dy- 
namics of liquid HF by inelastic x-ray scattering (IXS) 
at fixed temperature demonstrating the presence of both 
structural and microscopic relaxation processes. In the 
present paper we present an extended study of HF as a 
function of the temperature in the liquid phase between 
Tb = 292 K the boiling point and Tm = 193 K the 
melting point. Comparing the results obtained with two 
different techniques, IXS and Brillouin light scattering 
(BLS), we find the presence of a structural relaxation 
process in the entire explored temperature range. The 
relaxation time TaiT), in the sub-picosecond time scale, 
follows an Arrhenius temperature dependence with an 



activation energy strictly related to the number of hy- 
drogen bonds. The paper is organized as follow: Sec. II 
is devoted to the description of the experimental aspects 
related to the IXS and BLS measurements of the dynamic 
structure factor of HF. Sec. Ill is dedicated to the data 
analysis, in Sec. IV the main results are discussed and 
finally in Sec. V the outcomes of this study are summa- 
rized. 



II. SAMPLE ENVIRONMENT AND 
EXPERIMENTAL SET-UP 

High purity (99.9%) hydrogen fluoride has been pur- 
chased by Air Products and distilled in the scattering cell 
without further purification. The sample cell was made 
out of a stainless steel block, this material is well suited 
to resist to the chemical reactivity of HF. To allow the 
passage of the incident and scattered beam, two sapphire 
windows of 250fim thickness and 6 mm diameter, have 
been glued on two holder plates which have then been 
screwed to the body of the cell. An o-ring of parofluor 
has been applied between the window holders and the 
cell to guarantee a good tightness. The whole cell has 
been thermoregulated by means of a liquid flux cryostat 
DC50-K75 Haake. Further details of the sample cell will 
be described elsewher»i& 



A. Brillouin light scattering experiment 
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FIG. 1: Brillouin light scattering spectra of hydrogen fluoride 
at the indicated temperatures. 



The dynamic structure factor of HF in the GHz range 
has been measured by Brillouin light scattering using a 
Sandercock-type multi-pass tandem Fabry-Perot interfer- 
ometer characterized by high contrast (> 5-10^*'), resolu- 
tion (FWHM « 0.1 GHz) and a finesse of about 100. The 
wavelength of the incident radiation was A = 514. 5nm 
and the light scattered by the sample was collected in the 
back-scattering geometry {9 = 180°). The free spectral 
range (FSR) was set to 10 GHz, the integration time was 
approximately 2.5 s/ channel. The polarization of the in- 
cident light was vertical while the light scattered by the 
sample was collected in the unpolarized configuration. 
The aim of the present measurement is to determine the 
frequency position and width of the Brillouin doublets as- 
sociated to the propagation of the sound modes. As the 
relaxation time for HF ,in the investigated temperature 
range, is in the sub-picosecond region, we do not expect 
any evidence of the mentioned relaxation processes in 
the GHz range. Thus from the measured Brillouin peak 
position and width, it is possible to extract information 
about the adiabatic sound velocity cq and the kinematic 
longitudinal viscosity v^. 



B. Inelastic x-rays scattering experiment 

The inelastic x-rays experiment has been carried out 
at the very high energy resolution IXS beam-line ID16 
at the European Synchrotron Radiation Facility. The 
instrument consists of a back-scattering monochromator 
and five independent analyzers operating at the (11 11 
11) Si Bragg reflection. They are held one next to the 
other with a constant angular offset on a 6.5 m long ana- 
lyzer arm. The used configurationii, gives an instrumen- 
tal energy resolution of 1.6 meV full width half maximum 
(FWHM) and a Q offset of 3 nm~^ between two neighbor 
analyzers. The momentum transfer, Q, is selected by ro- 
tating the analyzer arm. The spectra at constant Q and 
as a function of energy were measured with a Q resolution 
of 0.4 nm^^ FWHM. The energy scans were performed 
varying the back-scattering monochromator temperature 
with respect to that of the analyzer crystals. Further de- 
tails on the beam-line are reported elsewhereiS,. Each 
scan took about 180 min and each spectrum at fixed Q 
was obtained by summing up to 3 or 6 scans. 
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FIG. 2: Stokes part of the Brillouin light scattering spectra 
of HF at the indicated temperatures. The fit (solid line) are 
superimposed to the data (open symbols). The curves are 
shifted on the y-axis one respect to the other. 
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FIG. 3: (a) Sound velocity co and (b) Kinematic longitudinal 
viscosity vl from Table |I] as calculated from the position and 
the width of the Brillouin light scattering excitations in back- 
scattering geometry (Q — 180'') (full circles). The straight 
lines represent the linear fit to the data. 



III. DATA REDUCTION 

A. Brillouin light scattering 

Unpolarized Brillouin spectra collected in a temper- 
ature range between 234 K and 283 K are shown in 
Fig. n The quantities of interest are the position and 
the width of the Brillouin peaks directly related to the 
sound velocity cq and to the kinematic longitudinal vis- 
cosity Vl of HF. In order to extract these two parameters 
the experimental data have been fitted in a limited region 
around the inelastic peaks with the function obtained by 
the convolution of the instrumental resolution K(lo) with 
a damped harmonic oscillator (DHO) function: 



(0(Q)2-c^2)2 + (2^r(Q))2 



(1) 



where = is the bare oscillation frequency and 

2r is approximately the full width at half maximum 



(FWHM) of the sound excitations. The results of the 
fitting procedure are reported in Fig. |2] superimposed to 
the experimental data. 

By exploiting the relations ^{Q) — cqQ and vl = 2r/27r, 
the adiabatic sound velocity cq and the kinematic longi- 
tudinal viscosity are obtained. The exchanged mo- 
mentum Q values are determined via the relation Q = 
^Y^sm(|), where n is the refractive index and 9 is 180° 
in the used scattering geometry. The temperature depen- 
dent refractive index, n{T), has been obtained by using 
the Clausius-Mossotti relation: 



1 



liTf 



-■KPn{T)a 



(2) 



where p„ is the number density and a the optical po- 
larizability of the HF molecule. The latter quantity, a, 
has been obtained using the values of n(T w 293/^)^^ 
and p{T « and assuming no temperature de- 

pendence for this parameter which turns out to be a = 
(9.606 ± 0.008) X^. The data for n(T) have then been 
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FIG. 4; IXS spectra of HF at fixed temperature in tlie low-Q region. Tlie raw data (lines + symbol) are plotted together with 
the corresponding experimental resolutions (dashed lines). 



TABLE I: Values of the frequency position z/q and of the width 
2T of the Brillouin doublet as obtained by the fitting proce- 
dure described in the text. The Q values, are also reported 
together with the sound velocities Co and the kinematic lon- 
gitudinal viscosities i/^. 



T 


vo 


2r/27r 





Co 




(K) 


(GHz) 


{GHz) 




(m/s) 


{cm? /s) 


234 


2.74 


0.24 


0.029 


600 


0.0058 


243 


2.51 


0.20 


0.029 


550 


0.0049 


253 


2.39 


0.14 


0.029 


530 


0.0034 


265 


2.27 


0.10 


0.029 


500 


0.0026 


269 


2.21 


0.09 


0.028 


490 


0.0023 


273 


2.15 


0.07 


0.028 


480 


0.0017 


278 


2.05 


0.05 


0.028 


450 


0.0011 


283 


1.97 


0.06 


0.028 


440 


0.0014 



obtained at each temperature by using a and p{T) whose 
expression is given by — 



velocity cq and the kinematic longitudinal viscosity 
are reported in Tab. I and shown in Fig.|21 The quantity 
co{T) follows a linear behavior characterized by a tem- 
perature dependence well represented by the equation: 



co{T) = cso + B ■ T 



(4) 



with Co = (1290 ± 40) m/s and B = (-3.0 ± 0.1) m/sK. 
The same procedure has been applied to derive the kine- 
matic longitudinal viscosity for which the linear fit 
provides a temperature behavior described by the rela- 
tion: 



(5) 



with C = (-13.5 ± 0.7) ctoVs and D = (1960 ± 
170) cm^K/s. All the values of Q, of the fit parameters 
and of the calculated i/^ and Cq, are reported in Tabled 



B. Inelastic x-rays scattering 



p{T) =po+A-T (3) 

with po = (1.616 ± 0.003) g/cm^ and A = -(2.25 ± 
0.01) • lO^'^ g/crrv'K. The derived values of the sound 



The IXS measurements, performed to probe the dy- 
namics of HF in the mesoscopic regime, are compared to 
the BLS results of the previous section in order to char- 
acterize the transition from the hydrodynamic regime to 
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FIG. 5: Dispersion curves at the indicated temperatures. The 
upper full lines are the linear fits to the high-Q data. The 
lower dashed lines indicate the adiabatic sound velocity as 
measured by Brillouin light scattering as shown in Sec. Ill AI 
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the mesoscopic. In this case the S{Q, lo) has been studied 
at four temperatures in the range 214-283 K at T = 214 
K, T = 239 K, T = 254 K and T = 283 K as a function of 
the wave vector Q. It has been varied between l-^15nm^^ 
for all the studied temperatures excepted for T = 239 K, 



where it has been selected between 2 31r 



Each 



energy scan took 180 min and each spectrum at fixed Q 
was obtained by summing up to 6 or 3 scans. We re- 
port in Fig. 0] an example of the measured spectra at the 
investigated temperatures and at the indicated momen- 
tum transfer (dotted line); they are compared with the 
instrumental resolution aligned and scaled to the central 
peak (full line). 



1. Markovian approach 

A first raw analysis of the spectra has been done in 
terms of the Markovian approximation in the memory 
function approachSi. In this approach the S{Q,uj) is ex- 
pressed as: 



SiQ,u;) = IiQ)- 



uJoiQ)M'iQ,Lo) 



~ LUoiQ)^ - ujM"{Q, cj)]2 + [uM'iQ, w)]2 

(6) 

where ujq{Q)'^ = {K bT / mS {Q))Q^ is the normalized sec- 
ond frequency moment of S{Q, uj), Kb is the Boltzmann 
constant, m is the mass of the molecule and M'{Q,uj), 
M"{Q,uj) are respectively the real and the imaginary 
part of the Laplace transform of the memory function 
M{Q,t). In the Markovian approximation, the decay of 
the memory function is faster than any system time scale 
and is modelled with a S function in the time domainal 
in such a way that Eq. reduces to a DHO. To fit our 
data we used this function plus a Lorentzian to take into 
account the finite width of the quasi elastic central peak 
. The detailed balance and the convolution with the in- 
strumental resolution have also been taken into account 
during the fitting procedure. One of the parameters we 
are interested in, is ^1{Q) which corresponds to the fre- 
quency of the sound modes. Its dispersion curve (i.e. 
its Q-dependence)is shown in Fig. [S] at low Q for the 
four analysed temperatures. The data show a common 
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FIG. 7: Q-dependence of the relaxation times Ta{Q) from 
the viscoelastic analysis at the indicated temperatures. A 
constant fit in the low Q region (dashed line) is also reported. 



FIG. 8: Q-dependence of the parameter r^(Q) (full circles) 
in liquid HF at the analyzed temperatures from a viscoelastic 
analysis. The full lines are the parabolic fits to the low-Q 
data. 



behaviour in the entire investigated T range, namely a 
linear dependence in the Q range 4 — 5^7 nm~^, with a 
slope corresponding to sound velocities higher than the 
adiabatic values Co measured by BLS and reported in 
the previous section. In addition for lower Q, the ap- 
parent sound velocity c{Q) determined by IXS seems to 
show a transition from the low frequency value Co to the 
higher value. This result provides the necessary infor- 
mation to extend what recently observed in liquid HF at 
T = 239 K^^ to a wider temperature range. The increase 
of c{Q) with increasing Q, in fact, is interpreted as due 
to the presence of a relaxation process, the structural re- 
laxation, already observed in HF at T = 239 K and still 
present in the whole explored temperature region. 



2. Viscoelastic approach 

The existence of a relaxation process with a charac- 
teristic time T in the range of the probed sound waves 
(i.e. such that fir ~ 1), as evident from the dispersion of 
i7((3), calls for a more refined choice of the memory func- 
tion with respect to the Markovian approximation. To 
describe the effect of this relaxation in S{Q,uj), we use 



the memory function based on the viscoelastic model. In 
this approach we describe a two relaxation process sce- 
nario with a memory function M{Q,t) given by the sum 
of an exponential decay contribution and a 5-functioni^ 



M(Q, t) = (Q)e-*/-° + V^{Q)m (7) 

where (Q) = [caoa(Qf - co{Qf]Q'^, is the strength 
of the a process and T^{Q) — A^t^((5). As in pure 
HF at r = 239 JsT— in fact, one expects a structural pro- 
cess described by an exponential decay and a microscopic 
process, very fast respect to the investigate time scale^S 
described by a (5— function. This approach has been suc- 
cessfully applied in the past to describe the dynamics of 
simple liquids^^ and liquid metals^''^''*"'^-^^. The thermal 
contribution in the memory function has been neglected 
being the value of the specific heats ratio 7 close to 1. The 
experimental data have been fitted to the convolution of 
the experimental resolution function with the dynamic 
structure factor model given by the combination of Eq. |B1 
and0 

The Q dependence of the fit parameters Cooa(Q), cq{Q), 
Ta{Q) and r^((3) is described in the following . 
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FIG. 9: (a) Behavior of the sound velocities in HF as a 
function of the temperature : co (open squares), Cooa (open 
circles), (b) Sound velocities ratio Cooq/cq (symbols); the 
straight line represents a constant fit to the data. 
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FIG. 10: T dependence of the low-Q extrapolation of the re- 
laxation times Ta{Q) of Figure[7|together with the Arrhenius 
fit (fuU line) of Eq.lSl 



In Fig.|Slthe Q behavior of the sound velocities is shown 
in the low Q region and for all the investigated tempera- 
tures. Both he infinite Cooa{Q) and zero co{Q) frequency 
limiting values, as obtained from the viscoelastic fit, are 
reported (open symbols). The comparison with the ap- 
parent sound velocity c{Q) — Vl{Q)/Q (full symbols) as 
derived from the dispersion curve of Fig. [S] shows the 
transition of c(Q) from cq{Q) to Cooa[Q)- The consis- 
tency between the two independent analysis strongly sug- 
gests that this transition is governed by the a-relaxation 
process in the entire temperature range. 

The Q dependence of the relaxation time Ta{Q) is re- 
ported in Fig. \7\ at the four investigated temperatures 
and in the Q range 1 15 nm~^. It shows a constant 
behavior, within the error bars, in the low Q region and 
a very week decrease at increasing Q as already observed 
in water^^ and many other systems-i. The fit of the data 
in the 1-^7 nm~^ range yields values Ta{Q 0) reported 
in Tab. In Fig. |S1 the Q dependence of the strength 
of the microscopic relaxation, r^((5), is reported at the 
four analyzed temperatures. The data show a quadratic 
behavior and have been fitted with a parabolic function 

The values of the parameter D are reported in Tab ITTl as 
a function of the temperature. 



TABLE II: Low Q values of the parameters used to calculate 
the kinematic longitudinal viscosity i^l(O) of Eq.|5| co(0) from 
Eq.2] Cooa(O) from the linear fit of Fig.O r(0) and F derived 
from the viscoelastic analysis as described in the text. 



T 


co(0) 


Cooa(O) 


Tc(0) 


D 


ul{0) 


(K) 


(m/s) 


(m/s) 


(ps) 


(cm^/s) 


{cm? / s) 


214 


650 


1220 


0.47 


0.0030 


0.007 


239 


580 


1080 


0.32 


0.0022 


0.004 


254 


530 


1010 


0.21 


0.0026 


0.003 


283 


480 


980 


0.17 


0.0026 


0.003 



ferent parameters analysed in the previous paragraphs. 
These parameters fully characterize the collective dynam- 
ics of our system in the mesoscopic regime. The values of 
cqj Cooa and Caoa/co at the four investigated temperatures 
are reported in Fig. |51 As shown in Fig. Etb) the ratio 
Cooa/ Co is temperature independent in all the explored 
T-range and it is close to two as in the case of water^— . 
The temperature dependence of the structural relaxation 
time in the Q ^ limit has been deduced from Fig. 
Here the low Q part of Ta{Q) has been fitted using a 
constant function. The obtained values are reported in 
Fig.Elon a linear scale as a function of the temperature. 
In the explored temperature range, the Ta(T) behavior 
is well described by the Arrhenius law (full line): 



IV. DISCUSSIONS 



Ta(T) = roe^BT (8) 



This section is dedicated to the discussion of the tem- with an activation energy Ea = (1.9 ±0.2) kcal/mol and 
perature dependence of the low Q behavior of the dif- tq = (6 ± 2)10~^^s. The temperature dependence in the 
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FIG. 12: Arrhenius plot for the relaxation time as obtained 
from viscoelastic analysis of the dynamical structure factor 
for water (diamonds)-- and hydrogen fluoride (circles); The 
dashed lines indicate the best linear fit to the data and their 
slop give an activation energy of 3.8 kcal/mol for water and 
1.9 kcal/mol for HF. 



limit Q ^ of the last fit parameter D is reported in 
Tab. ^ It yields values which appear to be temperature 
independent being D = (0.170 ± 0.025) meV/nm''^ — 
(2.6 ± 0.4) • 10~^ err? js. This result is consistent with 
previous findings according to which the microscopic re- 
laxation is a temperature independent process'^^. Using 
the low Q values T^ifS) of Tab. ^together with the low Q 
extrapolations of the other parameters (see Tab. llT|) . it is 
possible to calculate the kinematic longitudinal viscosity 
from the relationSi; 

VL = T^mcLM-cim + ^ (9) 

where co is the adiabatic sound velocity measured by 
Brillouin light scattering as discussed in previous sec- 
tion. The derived values for vj^ are reported in Fig. 1111 
they are found to be consistent with the hydrodynamic 
data reported in Tab. . This equivalence gives fur- 
ther support to the validity of the employed viscoelastic 
model. In the same figure we also report two viscosity 
data determined by molecular dynamics (MD) simulation 
simulations. A recent MD study of the transport coeffi- 
cients (longitudinal and shear viscosity, thermal diffusiv- 
ity and conductivity) of hydrogen fluoride^* provides two 
values for the longitudinal viscosity 77^, one at T = 205-?^, 
7?l(T = IQhK) = 0.91 • IQ-^g/cTOs and the other at 
T = 279i^, ^L{T = 279/^:) = 0.38 • V^-'^gjcms. They are 
reported in Fig. after rescaling for the density of Eq.|3| 
according to the relation vl{T) = r]L{T)/p{T), they are 
quite consistent with the experimental data. 
In Fig. 1121 we report, on an Arrhenius plot, the com- 
parison between the relaxation times for hydrogen flu- 



oride and for water^^. The activation energy found in 
water, constant in the examined temperature range, was 
Ea = (3.8 ± 0.6) kcal/mol while the one for hydrogen 
fluoride is Ea = (1.9 ± 0.2) kcal/mol as previously dis- 
cussed. It is worthwhile to relate the values of the acti- 
vation energies to the different networks present in the 
two liquids. While hydrogen fluoride forms linear chains 
with one hydrogen bond on average for each molecule, the 
preferred arrangement of water is the three-dimensional 
tetraedric structure with two hydrogen bonds for each 
molecule. If we indicate with uhb-h^o and uhb-hf 
the number of hydrogen bonds for H2O and HF respec- 
tively and Ea-H20: Ea^HF the activation energies for 
the two liquids, we see that they satisfy the ratio: 

Eg-H-jO ^ Eg^HF ^-j^Q^ 
nHB-H20 nHB-HF 

In previous studies on wates^ the activation energy 
has been associated with that of the H-bond (w 
5kcal/mol)^. The result of Eg. 1101 strengthens the idea 
that the structural relaxation process involves the H- 
bond networks of the system and it seems also to suggest 
that in this case the activation energy of the process is 
related to the number of H-bonds to make and break and 
not only to the strength of each H-bond. 

V. CONCLUSIONS 

We have presented inelastic Brillouin light and inelastic 
x-rays scattering measurements of liquid hydrogen fluo- 
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ride, a prototype of the class of hydrogen bonded hq- 
uid systems, in a temperature range comprised between 
214 K and 283 K. We demonstrated that the collec- 
tive dynamics of liquid HF is characterized by a struc- 
tural relaxation process in the sub-picosecond time scale. 
In the explored temperature region this relaxation pro- 
cess affects the collective dynamics in a Q range between 
l-^7nm~^. An accurate analysis in terms of the viscoelas- 
tic model in the memory function approach allowed to ex- 
tract and determine the temperature dependence of the 
parameters describing the dynamics at microscopic level. 
The relaxation time related to the structural relaxation 
process follows an Arrhenius temperature behavior with 
an activation energy Ea which, compared with the value 



previously measured in liquid water, enables to establish 
a connection between Ea and the number of hydrogen 
bonds per molecule of the specific system. 
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